Magnetic Flux Tuning of Spin Chirality in Mott Insulators with Ring Exchanges 
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A manifestation of the many-body Aharonov-Bohm effect in the magnetic-flux-tuned Mott insu- 
lators with three-spin and four-spin ring exchanges, presents as an effective tool to manipulate the 
ground-state spin chirality, such as, tune the magnitude continuously, switch an abrupt jump, or 
even reverse its sign. Such a mechanism is demonstrated explicitly in both quasi-one-dimensional 
ladders and two-dimensional lattices with triangles as elementary plaquettes. 

PACS numbers: 71.10.Hf, 75.10.Jm, 72.80.Sk, 75.50.Ee 



Introduction. — Mott msulators (Mis), as the paradigm 
of strongly correlated materials, have been commonly 
considered to have only magnetic properties at low en- 
ergies due to their spin moments. However^it has been 
demonstrated recently by Bulaevskii et al. [l| that, due 
to a certain form of charge fluctuations, geometrically 
frustrated Mis may exhibit electric orbital currents ac- 
companying spin textures with chirality. This notion of 
spin chirality itself has already been an intriguing topic 
for decades in quantum magnetism, superconductivity 
and anomalous Hall effect 2]. And the effective spin 
Hamiltonian of Mis may contain a linear coupling of the 
spin chirality to an external magnetic field, as proposed 
by Motrunich [3] when studying the organic compound 
K-(ET)2Cu2(CN)3 in which possible spin liquids with 
spinon Fermi surfaces are of particular interest [1, 0, [B, 01 ■ 
At the heart of the theories by Bulaevskii et al. and 
Motrunich [H, it is the three-spin ring exchange 
(3SRE) in addition to the Heisenberg antiferromagnetic 
(AFM) two-spin coupling, which is a specific form of 
charge fluctuations in Mis [1, |^. This multiple-spin 
exchange concept, initiated by Thouless now ap- 

pears as essential in various strongly-correlated sys- 
tems: bcc solid '^He 
graphite 11 
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crystals 
wires [16 



solid '^He films adsorbed on 
Il3ll . two-dimensional (2D) electron Wigner 
15| , and zigzag Wigner crystals in quantum 
Several experiments have confirmed the pres- 
ence of four-spin ring exchange (4SRE) in cuprates IT], 
18|. And loading cold atoms into optical lattices opens 
another avenue to design ring exchanges [l3| • 

Despite the above encouraging advances, one still has 
no quantitative understanding that how can we effec- 
tively manipulate the spin chirality in Mis with ring ex- 
changes. It is both of interest and timely to address this 
problem, and here we conduct such a study of a frustrated 
spin- 1/2 system with 3SRE and 4SRE modulated by a 
magnetic flux. Employing exact diagonalization (ED) of 
finite systems, we consider both a (two-leg) triangular 
ladder and a 2D triangular lattice geometry with peri- 
odic boundary conditions (PBCs), which are the simplest 
systems on which both 3SRE and 4SRE are possible. 



Beyond the weak-magnetic-flux regime, we explore the 
large parameter space systematically, and demonstrate 
that: at specific combinations of exchange interactions, 
varying the magnetic flux strength enable us to tune con- 
tinuously the magnitude, switch an abrupt jump, or even 
change the sign, of the ground-state (GS) spin chirality. 

Model Hamiltonian. — With the nearest-neighbor (NN) 
Heisenberg AFM coupfing, the 3SRE and 4SRE terms 
modulated by a uniform magnetic flux, the spin- 1/2 
model Hamiltonian in a triangular ladder/lattice reads 

(ij) ijk£A 



i'jkl 



(1) 



ijkleO 



where is the spin operator on site i. Pijk which de- 
fined as Pi23 : |cri, 0-2, ca) jcra, (Ti, cr2), is the cyclic 
permutation of the three spins sitting on a triangular pla- 
quette, and satisfies P{^3 = P/23 = ^321- And similarly 
P1234 : |o'i,(J2,cr3,cr4) ~* |cr4, CTi , (72, cs) is the cyclic per- 
mutation of the four spins sitting on a rhombus consisting 
of two elementary triangles. The three sums in Eq. ^ 
run, respectively, over all NN bonds, elementary triangles 
and rhombi. Contrary to '^He systems in which the "^He 
atoms are neutral, and similar to electron Wigner crys- 
tals [iOl , a magnetic flux through the exchange path can 
change the nature of the ring exchanges in Mis, owing 
to the Aharonov-Bohm (AB) effect, (j) is the magnetic 
flux treading a triangular plaquette, in units of ^o/Stt 
{(po — hc/e is the flux quantum). We focus on the phys- 
ical parameter space with J, K3 > and K4 > 0, and 
vary the ratio J/K3 and K4/K3 with the setting = 1 
(as an energy unit) in the following calculations. 

For our spin-1/2 case, Pijk can also be written in terms 
of electron operators as P123 = {c\aCif3){clf^C2j){cl^C3a), 
and satisfles i(Pi23 - P321) = -4Si • S2 x S3 = -4xi23, 
where Xijk represents the local spin chirality 0]. In pre- 
vious studies where a magnetic flux is absent, the 3SRE 
term itself favors ferromagnetism (FM) and can be taken 
into account just by modifying J and allowing both J > 
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FIG. 1: (color online). Intensity plots of the GS spin chirality 
Xo in the (j)-J parameter space (setting K3 = 1) of the 8x2 
triangular ladder at various fixed 4SRE strengths A'4's. 



and J < pj, DJ. 

e"^Pi23 + e-"^P32i = 

S3) — (-fl23 + -^321) 



However even for a small 0, since 

COS(/)(Pi23 + -P32l)-4sin0(Si-S2 X 

- 40X123, the magnetic flux couples 
linearly to the spin chirality in the low-0 limit P, 0, Q , 
and therefore could probably induce a non-zero chiral- 
ity density. Here we are concerned with the parameter 
space of (j) G [0;7r], and hence consider the many-body 
AB effect induced by a strong magnetic flux. 

Triangular ladders. — We firstly consider the ladder ge- 
ometry. We study the GS averaged (local) spin chirality 
Xo = —{Xijk)o (here a negative sign is added for conve- 
nience) by varying (f>, J and K^. The typical ED results 
for a triangular ladder of the size 8x2 are shown in Fig.[TJ 
In the absence of 4SRE {K4 = 0), xo is non- negative 
in the parameter region (j) G [0, tt], as shown in Fig. [Tfa). 
Note that xo has the symmetry xo{(t>) = ^Xo(^0) = 
— Xo(27r — (/)), which has been numerically confirmed. The 
whole (j)-J parameter space is roughly separated into two 
regions: the bottom left corner (J < 2.4, 0/27r < 0.17, 
and uniformly colored) with Xo = and saturated FM 
(S'tot = <5'max), and larger region with Xo > and spin- 
singlet GSs (S'tot = 0). The quantum critical line between 
these two regions has also been verified through tracking 
the non-analyticities in the GS energy function Eo{(j), J). 
Later, we will show that the quantum critical line does 
not depend appreciably on the ladder size. 

For the triangular ladders with PBCs, because of the 
identity P123+P321 = 2SrS2+2S2-S3 + 2S3-Si-|-l/2, the 
Hamiltonian with K4 = will reduce to H = —N cos<j) + 
Ji E (y'r^^ ■ ^1 + "^2 E (i^r -8^+4 sin fee A Si ■ X 
Sfc with Ji — J ~ A cos (j) and J2 — J — 2 cos 0, where the 
superscript "inter" ( "intra" ) corresponds to the effective 
interchain (intrachain) two-spin coupling Ji (J2). At the 
left boundary line of Fig. [^a) with = 0, there is a 



quantum critical point J « 2.4 corresponding to J2/ J1 = 
—0.25 which separates the saturated FM phase [21j and 
the dimer phase [l^]. And at the right boundary line of 
Fig. [ija) with = tt, the GS is also the dimer phase since 
J > and (j) = IT gives J2/J1 > 0.5 [23|. The features 
of these states will be displayed and discussed later by 
various correlation functions. 

In the presence of 4SREs, there are even more inter- 
esting behaviors of xo, as shown in Figs. [Hb)-(f) with 
five typical K^s respectively. At K4 = 0.2 [Fig. Wih)], 
the saturated FM region shrinks in the J direction while 
expand a little in the 4> direction. At K4 = 0.4 and 
K4 = 0.6 [Figs.[Hc) and (d)], at the center of saturated 
FM region, there appears a negative-xo region in which 
the GSs are spin singlets (S'tot — 0). When K4 is fur- 
ther increased to 0.8 [Fig. Wie)] and 1.0 [Fig. [T^f)], the 
negative-xo region continues to expand and occupies a 
significant portion in the (j)-J parameter space. 




FIG. 2: (color online). Triangular ladders: GS energy per 
site Eo/N (left) and GS spin chirality xo (right) versus (j), for 
various J's, 7^4 's and ladder sizes A'^'s. 

An intuitive analysis for K4 > is much more difficult, 
than that in the simpler case of K4 = 0. However, it 
should be noted that the 4SRE operators satisfy P1234 — 
P4321 = 5(Pi23 + P234 + P341 + P412 - H.c.) = 2i(xi23 + 
X234 + X341 + X412). Therefore, e'2*Pi234 + e-'^'^P432i - 

cos 20(Pi234 + P432I ) - 2 sin 20(xi23 + X234 + X341 + X412 ) ■ 

Due to the opposite signs and the different AB periods 
of 3SRE and 4SRE terms, the low-0-limit coupling coef- 
ficient and the portion of negative-xo region depend on 
the competitions between them. 
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In order to address the effects of ladder sizes, we com- 
pare three sizes of 8 x 2, 10 x 2 and 12 x 2. The mainly 
considered quantities are xo and the GS energy per site 
Eq/N. From Figs. [lja),(c),(e), we can see that the 
Eq {(t>)/N curves coincide well with each other for all three 
ladder sizes. And the xo('/') curves [Figs.mb),(d),(f)] also 
present us rather identical behaviors. All these results in- 
dicate that in the thermodynamic limit (iV — > oo), both 
Eo{(j))/N and Xo(</') will not deviate obviously from these 
finite-size results. 
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FIG. 3: (color online). Triangular lattices: Eq/N (left) and 
Xo (right) versus (f), for various J's, K^s, and ladder sizes A'^'s. 



defined by dij = (1 
spins as P12 : |o'i,(T2} — ^ |o'2,fi}), and this projector 
gives 1 on a singlet and on a triplet. The dimer-dimer 
CF between two bonds is D{r) = {dijdki) — {dij){dki), 
and Di{r) [D2{r)] for two parallel (non-parallel) rung 
bonds between two chains. The chiral-chiral CF between 
two triangles is defined as X{r) = (xijkXimn), and Xi{r) 
for two up-triangles (or equivalently two down-triangles), 
while X2{r) for an up-triangle and a down-triangle. Note 
that if the two triangles have some sites in common, X{r) 
may have a small imaginary part, and only the real part 
is plotted [H. 

At each case when ^tot = •S'max = N/2, xo always van- 
ishes, and meanwhile C(r) is surely positive and almost 
a constant at any range r. Now we would take a closer 
look at the three kinds of CFs of the other cases, and 
focus on the 12 x 2 triangular ladder as an example. 
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FIG. 4: (color online). Various GS correlation functions (see 
text) versus the range r of the 12 x 2 triangular ladders at 
various J's, ^'s and K4, — 0. 



Triangular lattices. — For the 2D triangular lattices, we 
focus on two cases with the sizes of 5 x 4 and 6x4 [Fig. [3] . 
In the absence of 4SRE terms {K4 — 0), similar to the 
previous case of ladders, the Hamiltonian of a 2D trian- 
gular lattice with PBCs will reduce to H = —N cos (f> + 
(J - 4cos(/)) 'Zi^j) Si ■ Sj + 4:Siii(j)J2ijke/\ ■ x ^k- 
From Fig. [21 we can see that Eo{(j))/N and Xo{4>) display 
quite similar behaviors resulting from the competitions 
between 3SRE and 4SRE terms as the ladders, such as 
the sign changes and abrupt jumps of xo- 

Long-range correlations and ordering. — We now turn 
to three kinds of correlation functions (CFs) in ladders. 
The first is the spin-spin CF, defined as C(r) = (Sj-Si+r), 
where r is the range (in units of the lattice constant) be- 
tween two sites along the chain direction and takes inte- 
ger (half-integer) values for intrachain (interchain) spin- 
spin CF CAr)_\C2ir)]. The other two CFs are defined 
as follows [71, ll3|. The dimer operator on a bond (i, j) is 



We first consider the simpler cases with only 3SREs 
(Fig. |4]with K4 = 0). For J = 1.0, tuning (j)/2TT across 
an FM quantum critical point at 0.125 [Fig. (UJa)], both 
Ci{r) and 6*2(7') exhibit that the GS consists of two- 
period FM domains with opposite magnetization, which 
is a remnant signature of long-range FM ordering in the 
FM region; tuning (/)/27r further to 0.22 [Fig. g^b)] at 
which Xoi^t') takes a maximum, Ci(r) [C2(?')] shows weak 
FM (AFM) correlations, and both Xi{r) and X2(r) re- 
veal nondecaying long-range correlations; when (jj is in- 
creased to TT [Fig. HI^c)], C(r)'s and X(r)'s show fast 
decaying behaviors, while Z3(r)'s reveal the long-range 
dimer ordering. For J = 3.0, = [Fig. [^d)], Ci(r) 
[6*2 (r)] shows strong (weak) AFM correlations because 
of J2 > Ji, and I?(r)'s also show slowly decaying cor- 
relations; tuning 0/27r to make Xoi4') take a maximum 
[Fig. He)] and then to 0.5 [Fig. |li;f)], there kinds of CFs 
resemble the J = 1.0 cases. 
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FIG. 5: (color online). The same as Fig. HI but with nonzero 
4SREs {Ki > 0). 



Next, we progress to much more complicated and inter- 
esting cases with 4SREs (Fig.IH). For (J^Ki) = (1.0,0.6), 
at = [Fig. EKa)], the nonzero K4 makes the long- 
range FM correlations destroyed, C(r)'s and Z?(r)'s ex- 
hibit nondecaying fluctuations although X(r)'s show fast 
decaying behaviors; tuning (/)/27r to 0.08 [Fig. Efb)], at 
which xoi4') takes a negative minimum, the C(r)'s ex- 
hibit FM domains and X{r)'s reveal long-range corre- 
lations; then tuning (f)/2Tr to 0.20 [Fig. (He)], at which 
Xoi(t>) jumps to a large positive value, X(r)'s show long- 
range correlations. For [J^K^) — (2.0, 1.0), tuning 0/27r 
to 0.11 [Fig. ini^e)] at which xo(0) takes a minimum, 
the X(r)'s show strong two-period fluctuations; tuning 
0/27r further to 0.45 [Fig. MS)] at which xo(</') drops 
steeply, Xi{r) displays strongly fluctuating correlations, 
and Ci (r) also shows strong intrachain AFM correlations. 
For {J,Ki) = (4.0,1.0), tuning cj) to 0.09 [Fig. E^g)] at 
which Xo(0) takes a local maximum, the X(r)'s show 
strong two-period fluctuations; when tuning (/)/27r to 0.47 
[Fig.ElJi)], we can see some behaviors resembling the pre- 
vious case in Fig. Elf). 

Summary and discussion. — For a spin-1/2 system in a 
triangular ladder/lattice with NN AFM coupHng, 3SRE 
and 4SRE, and a uniform magnetic flux </>, we can effec- 
tively manipulate the GS spin chirality xo, such as tune 
continuously the magnitude of xo by varying (j), switch 
an abrupt jump near an FM phase boundary, or even 
reverse the sign of xo^ and change the low-0-limit cou- 



pling coefficient. Various CFs discover the characteris- 
tic long-range correlations accompanying the tuned spin 
chiralities. Such a mechanism presents a peculiar mani- 
festation of the many-body AB effect on quasi-localized 
spins in Mis. This magnetic flux tuning of spin chiral- 
ity is expected be observed in 2D organic compound k- 
(ET)2Cu2(CN)3, quasi-lD and 2D Wigner crystals, and 
cold atoms in optical lattices with ring exchanges. 

This work was supported by NSFC of China (No. 
10904130) and State Key Program for Basic Researches 
of China (No. 2006CB921802). ED calculations were 
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